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Abstract. We monitored the BL Lac object S5 0716+714 
in the optical band during the period November 1994 - 
April 1995, which includes the time of a 7 -ray observa¬ 
tion by the Energetic Gamma Ray Experiment Telescope 
(EGRET) on February 14-28, 1995. The light curves in 
the R and B bands show fast fluctuations superimposed 
on longer timescale variations. The color index correlates 
with intensity during the rapid flares (the spectrum is flat¬ 
ter when the flux is higher), but it is rather insensitive to 
the long term trends. Over the 5 month observational pe¬ 
riod the light curve shows an overall brightening of about 
1 mag followed by a fast decline. The EGRET pointing 
covers part of the very bright phase (V ^ 13.2) and the 
initial decline. An ultraviolet spectrum was also obtained 
with the International Ultraviolet Explorer (lUE) (1200- 
3000 A) during the EGRET observations. The variabil¬ 
ity of the optical emission of S5 0716+714 by itself sets 
important constraints on the magnetic field strength and 
on the physical processes responsible for it. Interpreting 
the whole electromagnetic spectrum with synchrotron self 
Compton models leads to the prediction of a bright 7 -ray 
state during the EGRET pointing. We discuss how the 
7 -ray data could be used as a diagnostic of the proposed 
models. 

Key words: BL Lac objects: general - BL Lac objects: 
individual: S5 0716+714 - gamma-rays: theory - ultravi¬ 
olet: galaxies 
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1. Introduction 

One of the most important recent discoveries in the field 
of active galactic nuclei (AGN) is that flat spectrum ra¬ 
dio sources (blazars) emit a substantial fraction, some¬ 
times most of their power, at 7 -ray energies. About 60 
flat spectrum radio sources have been detected up to now 
above 100 MeV by the Energetic Gamma Ray Experiment 
Telescope (EGRET) onboard the Compton Gamma Ray 
Observatory (CGRO) (Fichtel 1994; von Montigny et al. 
1995; Thompson et al. 1995). Strong and rapid variabil¬ 
ity seems to be a common property of these sources in 
the 7 -ray band (e.g. von Montigny et al. 1995). The ori¬ 
gin of this high energy emission is still unclear. A point of 
agreement is that the emitting plasma is in relativistic mo¬ 
tion at small angles to the line of sight, with consequent 
beaming of the emitted radiation, which is required to 
avoid strong suppression of 7 -rays due to photon-photon 
absorption (McBreen, 1979; for application to blazars see, 
e.g., Maraschi, Ghisellini & Celotti 1992; Dondi & Ghi¬ 
sellini 1995). 

The correlated optical and 7 -ray variability observed 
in some objects (Wagner et al. 1995, Hartman et al. 1996) 
indicates that the same population of electrons could be 
responsible for the flux in both bands, and suggests that 
the optical emission is due to the synchrotron process, 
while the 7 -rays are produced by the inverse Compton 
process between highly relativistic electrons and soft seed 
photons. However, the origin of the seed photons, the lo- 
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cation and size of the emitting region(s), and the degree of 
relativistic beaming of the high energy radiation are un¬ 
known. According to one interpretation, the seed photons 
can be produced by synchrotron emission internally to the 
emitting blob or jet (Maraschi, Ghisellini & Celotti 1992). 
Another possibility is that the blob itself could illuminate 
a portion of the broad line clouds, whose reprocessed line 
radiation can dominate the radiation energy density in 
the blob (Ghisellini & Madau 1996). In another scenario 
the seed photons are produced externally, by the accre¬ 
tion disk (Dermer & Schlickeiser 1993), by the broad line 
region illuminated by the disk (Sikora, Begelman & Rees 
1994), and/or by some scattering material surrounding the 
jet (Blandford 1993, Blandford & Levinson 1995). Finally, 
a dusty torus could provide IR radiation to be scattered 
in the y-ray band, as suggested by Wagner et al. (1995). 

For plausible values of the parameters, the relativis¬ 
tic electrons which emit y-rays by the inverse Compton 
process emit synchrotron photons in the IR-optical-UV 
range. As a consequence, the two emissions are expected 
to be strongly correlated. Optical observations can there¬ 
fore be crucial to understand the nature of y-ray emission 
in blazars since they allow to derive information on both 
the seed photons and the radiating high energy electrons. 

Renewed interest in optical monitoring is motivated 
also by the discovery that in the radio domain the flux of 
some blazars varies by 10-30% on a timescale less than 
a day (Quirrenbach et al. 1992). The inferred brightness 
temperature, of the order of 10^® K, poses severe difficul¬ 
ties to theoretical models, since it cannot be reconciled 
with the Compton limit of 10^^ K for the commonly es¬ 
timated values of the beaming factor (see, e.g. Blandford, 
1990). Possible explanations are that these variations are 
due to regions with extremely high beaming or to coher¬ 
ent processes or that they are extrinsic, due to interstellar 
scintillation. Simultaneous optical and radio monitoring is 
thus critical to discriminate different models. 

In this paper we present the results of recent optical- 
UV observations of the BL Lac object S5 0716+714. This 
source is an ideal target for both the research fields out¬ 
lined above, since it has been detected several times in the 
y-ray band and has shown intraday variability at optical 
and radio wavelengths (Wagner et al. 1996). The opti¬ 
cal observations were carried out between Nov. 15, 1994 
and Apr. 30, 1995, the International Ultraviolet Explorer 
(lUE) pointings were on Eeb. 27 and March 1, 1995 and 
cover the EGRET-CGRO viewing period (1995 Eeb. 14- 
28). 

The outline of the paper is as follows: in Sect. 2 we 
review the available information on S5 0716+714 at all 
wavelengths; in Sects. 3 and 4 we describe and discuss the 
optical observations, taken from different sites and tele¬ 
scopes, and the UV observation obtained with lUE simul¬ 
taneously with the phase 4 EGRET pointing; in Sect. 5 
we model the overall emission spectrum with a homoge¬ 
neous synchrotron self Gompton (SSG) model and with a 


relativistic inhomogeneous jet model. In Sect. 6 we draw 
our conclusions. 

2. The BL Lacertae object S5 0716+714 

S5 0716+714 was discovered in the Bonn-NRAO radio 
survey (Kiihr et al. 1981) of flat spectrum radio sources 
with a 5 GHz flux greater than 1 Jy, and was identified as 
a BL Lacertae object, because of the featureless spectrum, 
by Biermann et al. (1981). From the stellar appearance on 
deep plates Stickel, Fried & Kiihr (1993) derived a lower 
limit of z > 0.2 on the redshift, while Schalinski et al. 
(1992) give z > 0.3. Proper motions of VLBI components 
have been observed (Schalinski et al. 1992), and would 
result in an apparent superluminal motion with Papp = 4.6 
if z = 0.3 and Hq = 50 km s“^ Mpc“^. 

A high dynamic range VLA 1.5 GHz radio map shows 
a double structure (Antonucci et al. 1986) resembling that 
of an FR H radiogalaxy seen end-on, but higher frequency 
(5 and 8.5 GHz) VLA maps do not show the double struc¬ 
ture (Wagner et al. 1996). The FR H morphology con¬ 
flicts with the suggestion that the parent population of 
BL Lac objects should be FR I radiogalaxies (see e.g. Urry 
& Padovani 1995 for a review) but it is consistent with a 
more general “unified” scheme whereby both FR I and FR 
H radiogalaxies have similar relativistic jets in their cores 
(Maraschi & Rovetti 1995). 

A further “anomaly” is that the observed ratio of the 
core to the extended radio flux, which is R = 4.7 at 5 
GHz and R ~ 1 at 1.5 GHz (Perley, Fomalont & Johnston 
1982), is at the low extreme of the R-distribution of the 
radio-selected BL Lac objects (see e.g. Ghisellini et al. 
1993 and references therein). 

The radio core has an inverted spectrum [a = —0.42 
between 1.5 and 5 GHz; Perley, Fomalont & Johnston 
1980; F{i^) (X The overall flat radio spectrum seems 
the result of the superposition of the core and more ex¬ 
tended and steeper components (see the decomposition of 
the radio spectrum in Eckart et al. 1986), even if the in¬ 
traday radio variability suggests that at least some of the 
radio emission may come from much smaller regions. 

The source was monitored in the mm band by Steppe 
et al. (1992, 1993) and Reich et al. (1993). At these fre¬ 
quencies (90-230 GHz) the flux varied by a factor 2-3 on 
timescales of the order of a month. As shown by Reich et 
al. (1993), the radio flux at all the monitored frequencies 
increased after Jan. 1992, the time of the first detection 
by EGRET-GGRO. 

In the far infrared it was detected by the Infrared As¬ 
tronomical Satellite (IRAS) at 12, 60, and 100 /rm (Impey 
& Neugebauer 1988). 

Despite S5 0716+714 is a very interesting source it 
has been very poorly observed in the past. Biermann et 
al. (1981) report the two photographic magnitudes 13.2 
and 15.5, the former derived from the POSS plate. The 
two previous visual magnitude estimates of 11 (Kiihr et al. 
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1981) and 19±1 (Perley, Fomalont & Johnston 1980) can¬ 
not be considered as fully reliable because no clear source 
for them is specified by these authors. The 8-month long 
relative R light curve, from September 1990 to March 1991 
(Heidt & Wagner 1996, Wagner et al. 1996), shows a varia¬ 
tion amplitude of about 2.5 mag. No historic curve, based 
on archive plates, to search for possible long term flux 
changes, is available. Takalo, Sillanpaa & Nilsson (1994) 
found a high level of optical polarization, which varied 
from 20% to ~ 16% in two days. 

In X-rays the source was observed by HEAOl, EIN¬ 
STEIN, and ROSAT (Biermannet al. 1981, 1992; Wagner, 
1992; Cappi et al. 1994; Comastri, Molendi & Ghisellini 

1995) . During the ROSAT observations the flux varied by 
70% in 1000 seconds and by a factor 7 in 2 days (Wag¬ 
ner 1992, Cappi et al. 1994). EGRET observed and de¬ 
tected S5 0716+714 several times from January 1992 (Lin 
et al. 1995; von Montigny et al. 1995). The integrated flux 
above 100 MeV varied between 1.3±0.5 and 5.3±1.3xl0“^ 
photons cm“^ s“^. Low statistics did not allow to derive 
a significant change of the energy spectral index during 
variations of the flux (Lin et al. 1995). Combining obser¬ 
vations of four different viewing periods, Lin et al. (1995) 
found an average energy spectral index a-f = 0.85 ± 0.20. 

Following the discovery of intraday radio variability 
with a ~ 10% amplitude (Heeschen et al. 1987), the 
source was monitored several times at radio wavelengths, 
always showing intraday variability (Quirrenbach et al. 
1989, 1992). At optical frequencies, variations on timescale 
of ~1 day with amplitude of several tenths of magnitude 
are reported by Wagner et al. (1996). 

A simultaneous optical {R band) and radio (5 GHz) 
monitoring campaign, performed for a month (February 
1990), suggested that the variations in the two bands were 
correlated, with amplitudes up to a factor 3 in the optical 
and 25% in the radio, and with typical timescales of 1 and 
7 days in both bands (Wagner et al. 1990; Quirrenbach et 
al. 1991). Further support to the suggestion of a correlated 
optical and radio variability comes from the correlation 
between the radio spectral index in the range 5-8.4 GHz 
and the optical flux (Wagner et al. 1996 and references 
therein). If the redshift 2 ; > 0.3 the implied brightness 
temperature exceeds 10^® K (Quirrenbach et al. 1991). 

3. Observations and data analysis 

3.1. Optical 

Photometry was carried out with several telescopes: the 
1.05 m astrometric telescope of the Observatory of Torino 
(B and R bands mainly), the 0.4 m automatic telescope 
of the University of Perugia (Tosti, Pascolini & Fiorucci 

1996) , the 0.5 m telescope of the Astronomical Station 
of Vallinfreda, and the 0.4 m telescope of the lAS-CNR 
at the Observatory of Monte Porzio (the last two oper¬ 
ated by the Rome group). All telescopes were equipped 


with GGD cameras and filters B, V (Johnson) and R, I 
(Gousins). Data reduction was performed with different 
software packages: the standard IRAF and MIDAS proce¬ 
dures and some others developed locally, as the ROBIN 
program developed at the Observatory of Torino by L. 
Lanteri (private communication). All images were bias and 
flat-field corrected. The comparison among data obtained 
in the same nights with different telescopes yielded the 
same results within the errors, which are typically 0.02 
mag in R and 0.03-0.05 mag in B. 

We did not find any published comparison sequence in 
the field of S5 0716+714; only one nearby star was cali¬ 
brated by Takalo, Sillanpaa & Nilsson (1994). We there¬ 
fore calibrated some reference stars using several Landolt’s 
(1992) fields. The stars are indicated in the map of Fig. 
1 and the corresponding magnitudes are listed in Table 1; 
quoted errors include the uncertainty on the absolute cali¬ 
bration, while magnitude differences between the reference 
stars in a given filter are good within 0.02 magnitudes. 
The star C of this sequence is the one already measured 
by Takalo, Sillanpaa & Nilsson (1994): our magnitudes are 
consistent with theirs. 


1 ’ • 



Fig. 1. Finding chart of S5 0716+714. See Table 1 for the 
magnitudes of the labeled reference stars. 


The reddening correction in the direction of S5 
0716+714 was taken Ay = 0.23. It was derived from 
the A^h value of 3.9 x 10^® cm“^ (Giliegi, Bassani & Ga- 
roli 1993), in agreement with the best fit values of both 
UV and X-ray observations (Cappi et al. 1994), using 
the conversion formula by Shull & van Steenberg (1985) 
and R=3.09 (Riecke & Lebofsky 1985). From the curve 
of Riecke & Lebofsky (1985) we derived also Ab = 0.30, 
Ab = 0.19, and Aj = 0.14. In the following we use these 
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values when deriving fluxes from the apparent magni¬ 
tudes. 


Table 1. Reference stars in the field of S5 0716+714 



I 

R 

V 

B 

A 

10.92+0.04 

11.21+0.04 

11.51+0.04 

12.07+0.06 

B 

11.79+0.05 

12.12+0.05 

12.48+0.05 

13.06+0.08 

G 

12.85+0.05 

13.18+0.05 

13.58+0.06 

14.17+0.08 

D 

12.97+0.04 

13.27+0.04 

13.66+0.04 

14.25+0.05 


The observation journal is reported in Table 2: here, 
for each date are given the measured magnitudes and their 
uncertainties in the four filters together with the used tele¬ 
scope. In some occasions several images were made in the 
course of the same night to search for intraday variability; 
these nights are indicated by an asterisk. The magnitude 
uncertainties have been evaluated by taking into account 
the mean dispersion of the internal error, but not the er¬ 
rors on the absolute calibration. The entire data set can be 
retrived via anonymous ftp at the address ftp.to.astro.it 
directory /pub/hlazars/. 

The light curves in the four bands are shown in Fig. 
2a,d; in Fig. 3 the part of the R light curve (the most sam¬ 
pled in that period) corresponding to the EGRET pointing 
is shown in more detail. 



650 700 750 800 850 

JD (2449.. + ) 


Fig. 2. Light curve of S5 0716+714 in the I, R, V, B bands. 


S5 0716 + 714 (14 Feb - 2 Mar 95) 

□ ^ TO <> ^ PG A ^ RM 



Julian Date - 2449000 

Fig. 3. Enlargement of the light curve in the R filter during 
the period of the EGRET pointing. 

3.2. Optical and Ultraviolet spectra 
— Optical spectrum — 

S5 0716+714 was observed with the 1.5 m Cassini tele¬ 
scope of the Loiano Station, operated by the Observa¬ 
tory of Bologna. The 4000-8000 A spectrum was obtained 
on March 19, 1995 at 00:58 UT, with the BFOSC all- 
transmitting spectrograph. A Thomson 1024 x 1024-pixel 
coated CCD detector was used with 300 1/mm grating and 
2 ”.5 slit to give a spectral resolution of ~ 25 A. The to¬ 
tal integration time was 2400 s. The 2-dimensional frame 
was biased and flat-field corrected using standard tech¬ 
niques, and for the extraction of the one-dimensional spec¬ 
trum the optimal extraction algorithm (Horne 1986) was 
adopted. The wavelength-calibration was performed with 
a fifth-order polynomial fitting to the dispersion curve 
of He-Ar comparison lamp spectrum. The spectrum has 
been flux-calibrated with the observation of the standard 
star BD +33 2642 (Stone 1977) and corrected with the lo¬ 
cal mean atmospheric extinction curve. The zero point of 
the spectral-flux calibration is not very accurate due to 
the poor seeing (~ 3 arcsec), and to the non-photometric 
conditions during the night. 

The resnlting spectrum is shown in Fig. 4, where no 
significant features are visible. The continuum spectrum 
was fitted with a power-law model assuming the redden¬ 
ing correction value Ay = 0.23 and the parametric formula 
for the mean extinction law computed by Cardelli et al. 
(1989), derived from the Riecke & Lebofsky (1985) extinc¬ 
tion curve. The best fit spectral index and associated Icr 
error is a = 0.95 ± 0.05. 
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Table 2. Observed magnitudes of S5 0716+714 


Date 

JD 

B 

V 

Rc 

Ic 

Tel. 

14 Nov 1994 

671.55 



13.52+0.03 


T 

15 

672.55 



13.73+0.02 


T 

23 

680.61 



13.86+0.02 


T 

29 

686.64 



13.65+0.02 


T 

11 Dec 

698.50 

13.87+ 0.03 




T 

02 Jan 1995 

720.51 

14.25+ 0.03 


13.41+0.02 


T 

03 

721.40 

14.15+ 0.03 


13.30+0.02 


T 

04 

722.52 

14.09+ 0.03 


13.26+0.02 


T 

06 

724.39 

13.84+ 0.03 


13.04+0.02 


T 

09 * 

727.39 

13.70+ 0.03 


12.86+0.02 


T 

14 * 

732.47 

13.96+ 0.03 




T 

24 

742.37 

13.80+ 0.03 


12.99+0.02 


T 

25 

743.33 

14.11+ 0.03 


13.25+0.02 


T 

30 

748.39 

13.77+ 0.03 


12.93+0.02 


T 

02 Feb 

751.44 

13.77+ 0.03 


12.93+0.02 


T 

04 * 

753.47 

13.90+ 0.03 


13.05+0.02 


T 

09 

758.45 

13.75+ 0.03 


12.90+0.02 


T 

14 

763.35 

13.73+ 0.03 


12.88+0.02 


T 

16 

765.38 

13.68+ 0.03 


12.85+0.02 


T 

18 

767.46 

13.64+ 0.03 

13.22+ 0.02 

12.78+0.02 


T 

19 

768.42 

13.63+ 0.03 

13.21+ 0.02 

12.78+0.02 


RT 

20 

769.25 

13.67+ 0.03 

13.23+ 0.03 

12.80+0.02 

12.29+ 0.03 

PR 

22 

771.37 

13.63+ 0.03 

13.25+ 0.04 

12.83+0.02 

12.31+ 0.04 

P 

27 * 

776.37 

13.88+ 0.03 

13.49+ 0.02 

13.07+0.02 


PT 

28 * 

777.44 

13.80+ 0.03 

13.42+ 0.02 

12.98+0.03 

12.45+ 0.03 

R 

01 Mar 

778.39 

13.94+ 0.03 


13.10+0.02 


T 

02 

779.37 

14.07+ 0.03 


13.24+0.02 


T 

05 

782.46 

14.29+ 0.03 

13.83+ 0.02 

13.45+0.02 

12.98+ 0.05 

R 

06 

783.30 

14.34+ 0.03 


13.51+0.02 


T 

07 

784.31 

14.30+ 0.03 

13.96+ 0.03 

13.48+0.02 

12.96+ 0.03 

P 

09 

786.41 

14.13+ 0.03 


13.32+0.02 


T 

10 

787.47 

14.17+ 0.03 


13.32+0.02 


T 

11 * 

788.39 

14.38+ 0.03 

13.98+ 0.03 

13.53+0.02 

13.01+ 0.03 

PT 

14 

791.31 

14.44+ 0.03 

14.00+ 0.03 

13.59+0.02 

13.07+ 0.03 

P 

17 

794.37 


14.41+ 0.03 

13.95+0.02 

13.40+ 0.03 

P 

19 

796.36 

14.81+ 0.03 

14.36+ 0.03 

13.97+0.03 

13.44+ 0.05 

R 

20 

797.26 

14.41+ 0.03 

14.08+ 0.02 

13.64+0.02 


P 


797.36 

14.34+ 0.03 

13.97+ 0.03 

13.54+0.02 

13.03+ 0.03 

R 

21 

798.30 

14.57+ 0.03 

14.15+ 0.03 

13.70+0.02 


T 

22 

799.37 

14.96+ 0.03 

14.52+ 0.02 

14.07+0.02 

13.52+ 0.03 

R 

23 

800.39 

15.06+ 0.03 

14.64+ 0.03 

14.20+0.02 

13.63+ 0.03 

P 

24 

801.37 


14.51+ 0.02 

14.05+0.02 

13.50+ 0.03 

P 

31 

808.39 

14.34+ 0.03 

13.97+ 0.03 

13.54+0.03 

12.98+ 0.04 

P 

01 Apr 

809.31 

14.47+ 0.03 

14.05+ 0.02 

13.62+0.02 

13.14+ 0.03 

R 

03 

811.40 

14.80+ 0.03 


13.85+0.02 


T 

04 

812.34 

14.71+ 0.03 


13.82+0.02 


T 

05 

813.31 

14.88+ 0.03 


13.96+0.02 


T 

06 

814.36 

14.77+ 0.03 


13.91+0.02 


T 

07 

815.33 


14.14+ 0.02 

13.70+0.02 

13.19+ 0.03 

P 

09 

817.33 

14.63+ 0.03 

14.25+ 0.02 

13.82+0.02 

13.28+ 0.03 

P 

12 

820.35 

14.94+ 0.03 


14.05+0.02 


T 

13 

821.45 

14.83+ 0.03 


13.98+0.02 


T 

14 

822.38 

14.76+ 0.03 


13.94+0.02 


T 

19 

827.33 


14.22+ 0.03 

13.80+0.02 

13.29+ 0.03 

P 

30 

838.38 

14.76+ 0.03 


13.95+0.02 


T 


The asterisk indicates the nights in which more than four observations with the same telescope were performed. 
Telescope identifications: R = Roma , P = Perugia, T = Torino 
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Fig. 4. The (binned) UV spectrum is plotted together with 
simultanoues optical data and with the optical spectrum made 
19 days after, as labelled. 

— Ultraviolet spectrum — 

S5 0716+714 was observed with the two cameras 
aboard lUE on February 27 (LWP) and on March 1, 1995 
(SWP, see Table 3). The spectral signal was extracted 
from the line-by-line image files using the GEX routine 
(Urry & Reichert 1988) running within MIDAS, and flux 
calibrated using the curves provided by Bohlin et al. (1990, 
SWP) and Cassatella et al. (1992, LWP); the resulting 
spectral flux distribution, where no apparent emission or 
absorption features are visible, is also given in Fig. 4. The 
flux was integrated and averaged over bins of 100 A, and 
the statistical uncertainties associated with the resulting 
flux points were evaluated as in Falomo et al. (1993). To 
the flux errors a 1.25% photometric error was added, ac¬ 
cording to Edelson et al. (1992). The SWP and LWP flux 
distributions were separately fitted to a power-law model 
taking into account the reddening correction Ay = 0.23, 
and the extinction curve of Cardelli et al. (1989). The 
fitted spectral parameters at the 1-a confidence level are 
given in Table 3. 

The optical photometry indicates that the source was 
stationary during the 2 days which separate the lUE ob¬ 
servations, thus, assuming absence of significant interday 
UV variability a power-law fit has been attempted to the 
joined SWP and LWP spectral flux distributions (1230- 
2800 A). The resulting spectral index is intermediate be¬ 
tween those obtained from the individual fits (a^ = 1.43 
± 0.03), with an acceptable value of Xr- This, together 
with the fact that the slope of the high energy part of the 
UV spectrum is not significantly larger than that of the 
lower energy region (Aa ~ 0.2), excludes the measurable 
presence of a curvature or of a bump in the spectrum at 

~2000 A. 

A power-law fit to the merged SWP and LWP spectra 
with the interstellar extinction taken as a free parameter 


yields a minimum for Ay coinciding with the galac¬ 
tic value. Including the simultaneous optical data (see the 
B, V, R and I data points in Fig. 4) the combined opti- 
cal+UV fit yields a =1.27± 0.02, with = 5.1 (16 d.o.f.). 
The BVRI data are fitted by a power law of slope a =1.1+ 
0.1, significantly flatter than the UV slope, therefore one 
can infer that there is a break in the spectrum which oc¬ 
curs around the U band. 

We finally note that the flux at 1400 A (J% = 3.21 ± 
0.07 mJy), is more than 20 times larger than that mea¬ 
sured in 1980 (lUE exposure SWP 9440; Pian & Treves 
1993). 

4. Results of the optical monitoring 

The light curves 

From the light curves of Fig. 2a,d, one can see that the 
source brightened after JD=2449720, remained bright for 
~50 days, and then declined rapidly, to reach a minimum 
around JD=2449800. Note that in the ‘high’ state the 
source almost reached its recorded historical maximum. 



740 760 780 800 820 840 


JD (2449..-h) 

Fig. 5. The light curve in the R band (panel a) and the adopted 
spline function passing through the local minima of the light 
curve (dashed line). In panel b) we plot the ratio (minus 1) 
between the observed flux and the spline value. The labels and 
the vertical dashed lines identify interesting flaring events. 

Superposed to this relatively long term behaviour, 
there are fast fluctuations, whose amplitude appears wider 
in the ‘low state’, indicating that they could contain the 
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Table 3. lUE observations 



AA 

MX) 

CXi/ 

Fin) (mJy) 

x7« 

d.o.f. 

SWP 54002 

1230-1930 

1600 

1.60+0.09 

3.98+0.07 

0.99 

6 

LWP 30122 

2200-2800 

2600 

1.39+0.23 

7.87+0.14 

0.25 

5 

SWP+LWP 

1230-2800 

2000 

1.43+0.03 

5.45+0.08 

1.01 

12 


AA indicates the wavelength range used in the analysis. 


1 (T errors are given both for the spectral index and for the flux. 

The monochromatic flux F{v) is calculated at the frequency corresponding to the wavelength A. 


same amount of energy as in the ‘high state’. To make this 
behaviour more apparent we modelled the ‘long’ trend by 
interpolating a smooth curve through the local minima of 
the R light curve. We used a cubic spline interpolation, 
which is plotted as dashed line in Fig. 5a. We stress that 
this spline interpolation crucially depends on the chosen 
minima of the light curve, and therefore it should be con¬ 
sidered only as a possible representation of the source be¬ 
haviour. The percent excesses between the observed fluxes 
and those derived from this interpolating curve are shown 
in Fig. 5b: labels from A to E correspond to five events of 
’short’ timescale. The event C is the best sampled, whereas 
the shape of the event B can depend on the adopted spline 
and could be considered as the superposition of two (or 
even more) events. Notice that the event A has the short¬ 
est duration and the smallest variability amplitude. For 
the best sampled peak in this curve, the rising and decay¬ 
ing times appear approximately equal, with no indications 
of a ‘plateau’ longer than one day. 


4-2. Intranight Variations 

Rapid variations of the flux of S5 0716+714 have been ob¬ 
served in different occasions (e.g. Wagner & Witzel, 1995). 
In our campaign we also searched for intranight variations 
especially in the R band (see Table 2). Variations with an 
amplitude of about 5% in few hours in the course of the 
same night were observed in some occasions. Some of the 
resulting light curves in the B and R bands are shown in 
Figs. 6-7. The lower panels in these figures show the dif¬ 
ference between the instrumental magnitudes of the two 
reference stars used in the night (normalized to the mean 
value). 

Notice that in the two nights of 1995 Jan. 9 (Fig. 6) 
and Feb. 4 (Fig. 7) evident trends in the source brightness 
are apparent, like if our sampling was limited to a fraction 
of a longer time scale change. These cases, however, are 
important because the probability of their occurrence can 
be evaluated by means of the run test - see, for instance, 
the introductory textbook by Barlow (1989) - which is 
non-parametric and independent of the measured uncer¬ 
tainties. In the night of Jan 9, 1995, we have 9 data points 
in the B band and other 10 in the R band, both show¬ 
ing decreasing brightness trends and having each one two 
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Fig. 6. Intranight variations during the Jan. 9 1995 night, in 
the R band (left panel) and in the B band (right panel) 
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0.40 0.42 0.44 0.46 0.48 0.50 0.52 0.54 
Julian Date - 2449753 
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Fig. 7. Intranight variations during the Feb 4 1995 night, in 
the R band (left panel) and in the B band (right panel) 
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runs only. The expected run numbers for these two data 
samples are 5.44 {B band) and 6.0 {R band), correspond¬ 
ing to 2.5 and 2.7 standard deviations, respectively. Their 
probabilities, estimated from a one-tailed gaussian distri¬ 
bution are 6.2 x 10“^ and 3.7 x 10“^; the combined chance 
probability to observe these trends at the same time and 
in the same direction in the two bands can be estimated as 
1/2 times their product, because the observations in the 
two bands can be considered independent, as indicated by 
the comparison star behaviour. We obtain then the value 
of 1.2 X 10“®, which gives a high statistical significance to 
this effect. 

Again on Feb 4-5 a well defined trend in the same 
direction is apparent in the two bands. The run test gives 
a total probability of 5.9 x lO”'^, greater than the previous 
value because of the smaller number of data points, but 
statistically significant. 

4-3. Flux-spectral index correlation 

We computed the spectral index a (by a fit procedure), af¬ 
ter dereddening the observed fluxes, for the nights in which 
at least three filter data are available, including intranight 
data. These nights cover the period from February 18 to 
the beginning of April in which the source brightness de¬ 
clined by 1.4 magnitudes. The resulting values of a vary 
between 0.81 and 1.15, with a mean value of 0.94, and 
do not have any significant correlation with the flux: the 
linear correlation coefficient with the V magnitude is only 
0.15. In particular the spectral slopes measured when the 
source was at the maximum (February 19) and minimum 
(March 23) brightness were both equal to 0.99. 

We also computed the 2-point spectral index aBR^ 
which can be compared with the spectral index calculated 
using more filters. The differences between these indices 
are usually less than the statistical uncertainties and this 
gives confidence that ubr traces the ‘true’ optical spectral 
index. Again these values vary in the range 0.83-1.2. 

We found that when the source reached the “low” state 
(after JD=2449790) the spectral index responded to the 
fast flux changes in the sense that the spectrum was flatter 
when the source was brighter. This is evident from Fig. 8 
where R and aBR are plotted against time. The linear cor¬ 
relation coefficients of aBR with B and R were 0.627 and 
0.477, respectively, and the corresponding chance proba¬ 
bilities are 7.1x10“^ and 5.3x10“^. As pointed out by 
Massaro & Trevese (1996) the latter estimate is biased to¬ 
wards a lower significance of the correlation and therefore 
this probability must considered as a conservative indica¬ 
tion of the correct significance. A less biased estimate of 
the true correlation coefficient can be obtained by taking 
the mean of the two previous values, and in this case the 
corresponding probability results equal to 2.1x10“^. We 
can conclude that S5 0716+714 shows a significant (> 95 
%) positive correlation between the spectral index and the 
flux only during rapidly varying changes. 



800 820 840 

JD (2449.. + ) 


Fig. 8. The light curve in the R band during the ‘low’ state 
of S5 0716+714 (panel a) compared with the corresponding 
spectral index obh (panel b). 


5. Discussion 

5.1. Short time variations 

The structure of the light curves of S5 0716+714 together 
with the color index behaviour suggest that two different 
processes are operating in the source: the first one is re¬ 
sponsible for the long term “achromatic” rise and decay of 
the brightness, while the second one is responsible for the 
fast flares for which the color index (aBn) is significantly 
bluer at higher intensity. It is possible to interpret the first 
behaviour in terms of emission produced in a large region, 
in which a stable energy injection mechanism operates, at 
least on few month timescale. 

The fact that the rising and decaying timescales of 
the short term flares are comparable, with much shorter 
plateaux, bears important consequences. In fact, it is 
highly unlikely that these two timescales correspond to the 
particle acceleration and cooling timescales, respectively, 
since a priori there is no reason for them to be equal, and 
the observed equality would be coincidental. Instead, the 
equality is naturally explained in two frames: 

1) As Wagner et al. (1993) suggested for the similar be¬ 
haviour of S4 0954+658, the characteristic variability pat¬ 
tern of the flares could be caused by a curved trajectory of 
the relativistic emitting blob; in this case the viewing an¬ 
gle 6 (i.e. the angle between the blob velocity vector and 
the line of sight) varies with time and consequently the 
beaming or Doppler factor <5 = [r(l — /3cos0)]“^ changes 
too. Since the observed monochoromatic flux can be as¬ 
sumed to be proportional to symmetric 9 variations 
can account for the observed flare shape. In this case, for 
rather small amplitude changes, one has 
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—■— — (3 + ct)—— (3 + Q!)/3r(5 A0 (1) 

r 0 

which shows that, for F ~ ~ 1 / siii 0 ~ 10 , a flux varia¬ 

tion of 20% - see the results of Sect. 5.2 - implies a change 
of 0 of less than 0.5 degrees only. 

2) The rising and decay times are both associated with 
the light crossing time. This can happen if both the elec¬ 
tron injection and cooling processes operate on time-scales 
shorter than R/c (where R is the source radius), since in 
this case the light crossing time is the only characteristic 
time that can be observed. On the other hand, the short 
duration of a plateau also requires that either the accel¬ 
eration or the cooling times (or both) are not much less 
than R/c. Assume in fact that both these times are much 
shorter than R/c. Then at any given time the received 
flux originates in a different partial volume of the source. 
Different times correspond to different sub-volumes of the 
source and, if these regions are equal, the observed flare 
would have a plateau. 

The situation is further complicated by the contribu¬ 
tion of the steadier component, which dilutes the flux and 
spectral index variations. 

We can envisage a test to discriminate between the 
two possibilities. In case 1) the variable component should 
have a constant spectral index: amplification due to dif¬ 
ferent degree of beaming would in fact be ‘achromatic’ (if 
there is no spectral break within the considered bands). 
The sum of the variable (flatter) and the steady (steeper) 
components results in a continuous flattening of the spec¬ 
trum as the source brightens. The flattest spectral index 
should correspond to the maximum flux. In case 2) the 
variable component has instead a variable spectral in¬ 
dex: the maximum flux corresponds to integration over 
the source volume, whose sub-regions contain electrons 
of different ages. Therefore the maximum flux does not 
correspond to the flattest spectral index. 

Our data are not sufficient to discriminate between the 
two possibilities, since a better sampling (in more bands, 
including IR and, possibly, UV) is needed. Nevertheless, if 
possibility 2 ) turns out to be preferred by future investiga¬ 
tions, it will be possible to draw important consequences 
about the value of the magnetic held strength and the par¬ 
ticle energy. We will derive these two parameters in the 
next section by independent arguments, and it is interest¬ 
ing to see if they agree. 

Assume therefore that the rapid flares are due to 
quasi-istantaneous injection of electrons and rapid cool¬ 
ing. In particular we require that the synchrotron and self 
Compton cooling times are equal to or shorter than R/c. 


^cool 


Girmc^ 


arcjoB^i^ + Us/Ub) 


< 


R 


( 2 ) 


where ax is the Thomson scattering cross section, 'yoiric^ is 
the energy of the particles emitting in the optical, Ub = 
B^/{8Tr) is the magnetic energy density, and Us is the 


radiation energy density of the synchrotron emission (as 
measured in the comoving frame). The observed frequency 
is given by 

Vo = -"foVB——, (3) 

3 1 + z 

where vb = eB/{2'Kmc) is the cyclotron frequency. By 
relating the radius R of the source with the observed vari¬ 
ability timescale through R ~ ctvar5/{l + z), we obtain 
the two limits: 




1/3 


7o < 2 X 


td(l + Us/U b){^ + z) 


1/3 


( 4 ) 

( 5 ) 


where B is in Gauss, td is the variability timescale mea¬ 
sured in days, and Vo = Hz. Note the relatively 

weak dependence on the beaming factor 5. Assuming 
(5/(1 + z) = 10 , td = I day and Us/Ub = 1 (as indicated 
by the ratio between the self-Compton and synchrotron 
luminosities), we find B > 0.2 G and 7 o < lO'^'. 


5.2. Homogeneous Synchrotron Self Compton Models 

In Fig. 9 we show the overall spectrum of S5 0716+714 
from radio to 7 -rays, calculated assuming a redshift z = 
0.3. Solid symbols in the optical {B, V, and R) and UV 
bands correspond to observations simultaneous (Feb. 27 
and Mar. 1, 1995) with those of the CGRO pointing, whose 
data are not yet public. Empty squares are not simulta¬ 
neous data taken from the literature. Empty circles refer 
to simultaneous radio to UV data published by Wagner et 
al. (1996), the X-ray ROSAT spectrum is from Cappi et 
al. (1994) and the 7 -ray data are from Lin et al. (1995). 

The spectrum of S5 0716+714, like that of other 7 - 
ray loud blazars, is characterized by two broad peaks in 
iy-iyF{iy) plots, which can be interpreted as due to the 
maximum output of the synchrotron and Compton com¬ 
ponents. Note that the knowledge of both the optical and 
UV spectra is crucial to define the frequency of the syn¬ 
chrotron peak. 

Among the various “Compton” models (see Introduc¬ 
tion) we focus on synchrotron self Compton (SSC) mod¬ 
els, since no emission line has been observed in this source. 
This limits the importance of photons produced externally 
to the jet or emitting blob, which could serve as seed pho¬ 
tons to be Comptonized in the 7 -ray band. 

Consider a spherical blob, having the variability 
timescale radius R = ctvarS/{l + z) (in the comov¬ 
ing frame), with a tangled magnetic field B, emitting 
synchrotron and self Compton radiation. For appropri¬ 
ate particle distribution (to be better specified later), 
the synchrotron and self Compton spectra peak at the 
observed frequencies bs (~ Bq), given by Eq. (3), and 
Be = ( 473 ) 7 /^ 5 . From the ratio bq/bs we derive 
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Fig. 9. Overall spectrum of S5 0716+714. Empty circles are 
from Wagner et al. (1996), ROSAT X-ray data are from Cappi 
et al. (1994) and 7 -ray data are from Lin et al. (1995). Filled 
squares are optical and UV data from our observations, simul¬ 
taneous (Feb 27 - Mar 1, 1995) with EGRET observations (not 
public yet). Empty squares are not simultaneous data from 
the literature. Lines refer to homogeneous SSG models, as dis¬ 
cussed in the text. 


B5 ~ 


v 


2 

s 


2.8 X 


(1 + ^) 


( 6 ) 


The ratio of the luminosities of the two peaks mea¬ 
sures the relative importance of the synchrotron and self- 
Compton powers. Considering the Compton scattering 
process occurring entirely in the Thomson regime, this 
ratio is directly related to the ratio between the magnetic 
and radiation energy densities of the synchrotron photons: 


Then 


.. Ub 

AttR^cUbS'^ 

(7) 

~ Us 

Ls 

^ (l + ^)( 

' 2L| +'> 

(8) 

^LccHl^J 


Equations (6) and (8) can be solved, yielding 

/ X 1/2 / r2 \ 1/'* 

(5 ~ 1.67X 10^ 


B 


ttiFs J \ -^C,45 

3 , 1/2 / y \ 
( LcA5 


-2.14x10-^\1 + z)^ ^ 


'c 


S,45 


(9) 

( 10 ) 


where Ts ,45 and Tc ,45 are measured in units of 10"^® erg 
s“^ and td is tyar in days. Despite their simplicity, Eqs. 
(9) and (10) are a very powerful tool for deriving B and 
S, but we must stress that for a reliable estimate we need 
to determine with great accuracy the simultaneous values 
of VC and vs, together with an estimate of the variability 
timescales. 

For S5 0716+714, in the high state we have Ls = 2.2 x 
10^^ erg s“^, Lc = 4.4 x 10^^ erg s“^, vs = 10^® Hz, 


VC = 8 X 10^^ Hz, and tyar = 2 days. With these values 
we derive B ~ 0.5 Gauss, S ~ 11, and jo — 7.7 x 10^. 
Notice that these B and 70 values are in agreement with 
the limits obtained in Sect. 5.1. 

To calculate the emission spectrum we need to know 
the energy distribution of the emitting electrons. In a sim¬ 
ple model with continuous injection of Qij) oc 7 “® [cm“^ 
s“^] relativistic electrons throughout the source, the rela¬ 
tivistic electron density N{j) can be obtained by solving 
the well known steady continuity equation, leaving as free 
parameters the injected power Linj, the spectral index s 
of the injected distribution and the extremes of the energy 
range 7 mm and jmax of the injected electrons (see Ghis¬ 
ellini 1989 for details). For s = 2.6, the resulting 1 V( 7 ) can 
be approximated with a double power law, with indices 2 
and 3.6 below and above 7 mm, respectively. Therefore we 
can identify 7 mm as 70 . The emission spectra, computed 
for low (dashed line) and high (solid and dotted lines) lu¬ 
minosity states are shown in Fig. 9. The two states corre¬ 
spond to an increase of the injected power from 10 ^^ 
erg s“^ to 3.3 x 10"^^ erg s“^. The values of the other model 
parameters were R = 3x 10^® cm, <5 = 12, ^min = 3 x 10® 
(or 4 X 10®, dotted line), 'ymax = 10®, and B = 0.48 (or 
0.84, solid line) G. The two high state spectra differ for 
the values of B and 7 mm: the one represented by the solid 
line has an increased magnetic field to maintain the ratio 
Us/Ub constant and close to unity, while in the other the 
magnetic field is assumed to remain constant, and conse¬ 
quently the ratio Us/Ub increases. 

These two choices, both reasonable from a physical 
point of view, yield 7 -ray luminosities which are a fac¬ 
tor of 3 to 10 greater than the EGRET luminosities of 
1992-1993. Furthermore, the emission is expected to vary 
simultaneously at all frequencies from the mm band up to 
the GeV band, although the amplitude of the synchrotron 
variations can be smaller than the Gompton one. 

5.3. Inhomogeneous jet model 

We adopt the model of Ghisellini & Maraschi (1989) (see 
also Ghisellini, Maraschi & Treves 1985; Maraschi, Ghis¬ 
ellini & Gelotti 1992) in order to check if relaxing the hy¬ 
pothesis of homogeneity changes substantially the inferred 
physical parameters. In this jet model the density of the 
emitting particles and the magnetic field strength decrease 
along the jet as power law functions of the distance x from 
the jet apex. Relativistic particle are assumed to have the 
same power law distribution, all along the jet, but the high 
energy cut-off (and the maximum emitted frequencies) are 
function of x. Spectral break and different spectral indices 
in the radiation spectrum are due to the different superpo¬ 
sitions of the local spectra. The inner part of the jet has 
a parabolic shape, with cross sectional radius r oc 
while the outer part is conical (r oc x). In addition, the 
emitting plasma flowing in the jet is assumed to accelerate 
outwards in the parabolic region, while it has a constant 
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bulk velocity in the conical part. We also fixed the depen¬ 
dence of the magnetic field with the cross sectional radius 
r of the jet, setting it to be B oc and the shape of 
the particle distribution oc 7 “^, resulting in a fixed 

local spectral index, Uo = 0.5). We furthermore required 
that the photon-photon opacity must not be important 
all along the jet and that the region producing most of 
the optical can vary in a ~1 day timescale. With all these 
constraints, the set of parameters which better interpolate 
the data is rather well defined. These are listed in Table 4, 
and correspond to the spectra shown in Fig. 10. We also 
plot, for the low state, the emission coming from different 
regions of the jet. 



Fig. 10. The overall spectrum of S5 0716+714 interpreted as 
emission from an inhomogeneous relativistic jet. The solid lines 
refer to the emission from the entire jet, in the ‘high’ and ‘low’ 
state, according to the parameters listed in Table 4. The syn¬ 
chrotron and the self Compton components are shown sepa¬ 
rately. Dashed and dotted lines refer to the emission from re¬ 
gions within X = SOOxo and within x = SOxo from the jet apex, 
respectively. See Table 4 for the definition of xo- 

The bulk of the 7 -ray emission (up to a few GeV) is 
produced where the parabolical and the conical part of the 
jet connect, at distance a; ~ 5 x 10^^ cm from the jet apex, 
where the cross sectional radius of the jet is r ~ 2 x 10 ^® 
cm. As explained in Ghisellini, Maraschi & Treves (1985) 
and in Ghisellini & Maraschi (1989), this is expected in 
this class of models, due to the dependence of the Comp¬ 
ton flux as a function of distance x, which is different in the 
parabolic and the cone zones. However, as can be seen also 
in Fig. 10, higher frequencies (both in the synchrotron and 
the self Compton branches) are produced closer in. This 
is not the consequence of the particular set of parameters, 
but is again a general characteristic of the model: requir¬ 
ing to fit the spectral change in the IR-UV-soft X-ray 
spectrum as a superposition of local spectra (each with 
a different high energy cut-off), inevitably implies higher 
frequencies be produced in the innermost regions of the 
jet. This characteristic can well be tested with future 7 - 


ray satellite, such as GLAST, able to observe up to the 
100 GeV band, where the emission is predicted to vary 
more dramatically and with the shortest timescales. Less 
dramatic differential variability is however expected also 
in the different EGRET energy bins. In the intermediate 
region between the paraboloid and the cone the magnetic 
field strength is H ^ 0.5-1 G, the bulk Lorentz factor is 
F cs 12, the corresponding beaming factor is <5 ~ 14, and 
the magnetic and radiation energy densities are close to 
equipartition. 

In conclusion, we derive a set of parameters which 
agree well with the findings of the homogeneous model. 
With two important differences: 

1) In the inhomogeneous jet model the synchrotron X- 
ray emission comes from the inner part of the paraboloid, 
at distances x ~ 10 ^^ — 10 ^® cm from the jet apex (and 
cross sectional radii r ~ 10^^ — 3.5 x 10^^ cm. Correspond¬ 
ingly, the X-rays should vary with the minimum timescale 
(the observed tyar ~ 1 hour), possibly explaining the very 
fast variability reported by Cappi et al. (1994). 

2) In the inhomogeneous jet model variability is si¬ 
multaneous only in selected bands. This can be seen in 
Fig. 10 where the dashed and dotted lines correspond to 
the synchrotron and the self-Compton emission arising in 
same regions of the jet: ~ 1 GeV emission should vary fast 
and simultaneously with 1-10 keV X-rays, while ~ 100 
MeV 7 -rays should vary simultaneously with the optical- 
UV emission and be characterized by a somewhat longer 
timescale. Spectral index variations in the IR-UV are pos¬ 
sible if a perturbations travelling along the jet activates 
different regions, starting from the innermost one, respon¬ 
sible for the X-rays (for a detailed discussion, see Celotti, 
Maraschi & Treves, 1991). In this case the higher frequen¬ 
cies vary first, followed by the lower ones. Note that the 
predicted correlation between the flux and the spectral in¬ 
dex is in the same sense as observed: flatter slopes when 
the flux is brighter. 

Finally, note that in the inhomogeneous jet model, as 
well as in the homogeneous one, one cannot explain simul¬ 
taneous intraday variability: even in the favorable case of 
a perturbance travelling down the jet, close to the line of 
sight, one should observe a time delay between the optical 
and the radio emission (even if shortened by relativistic 
effects). 

6. Conclusions 

Dense optical monitoring is a very powerful tool to better 
understand the radiation processes operating not only in 
the optical band, but also in the high energy part of the 
spectrum, and to put models for the 7 -ray emission on 
test. 

By themselves, the optical data we have obtained for 
S5 0716+714 are indicative of a complex behaviour possi¬ 
bly characterized by more than one component. The de¬ 
pendence of the spectral index on the fast variations of 
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Low state 
Paraboloid 

Low state 
Cone 

High state 
Paraboloid 

High state 
Cone 

Note 

Xo (cm) 

IQi® 

5 X 10^’’ 

10®® 

5 X 10®® 

initial distance 

Xmax (cm) 

5 X 10^’’ 

10i» 

5 X 10®'® 

10®® 

final distance 

Bq (Gauss) 

12 

0.54 

12 

0.54 

initial magnetic field 

No 

3 X 10® 

2.8 X 10® 

6.7 X 10® 

6.3 X 10® 

initial electron density 

ao 

0.5 

0.5 

0.5 

0.5 

local spectral index 

m 

1 

1 

1 

1 

B oc r-™ 

n 

1.5 

1.5 

1.5 

1.5 

N oc r-" 

a 

0.4 

0 

0.4 

0 

P oc 

'ymax 

4.2 X 10® 

1.4 X lO"' 

4.2 X 10® 

1.4 X 10“® 


Po 

1 

12 

1 

12 

initial bulk Lorentz factor 

6 (degree) 

4 

4 

4 

4 

viewing angle 


Table 4. Input parameters for the inhomogeneous jet model. Note: 2 = 0.3 is assumed. For the parabolic part, the cross 
sectional radius r scales with distance a: as r oc 


the flux and its insensitivity to longer-term trends, sug¬ 
gest localized injections of energy in the form of fresh 
and rapidly cooling relativistic particles, superimposed to 
a steadier emission, possibly coming from larger regions. 
Other constraints can be obtained by the shape of the 
fast flares, which seem symmetric and without plateau. 
In particular, this allows to derive a lower limit on the 
magnetic field strength. It is very interesting to verify if 
other sources behave like S5 0716+714 in this respect, to 
assess if this is a fundamental feature to be explained by 
emission models for blazars in general. Indeed, some of the 
best monitored surces (e.g. 3C 66A, Takalo et al. 1996; OJ 
287, Sillampaa et al. 1996; 0954+658, Wagner et al. 1993; 
PKS 2155-304, Pian et al. 1997; PKS 0422+004, Massaro 
et al. 1996), seem to confirm this behaviour. 

S5 0716+714 was very bright in the optical during the 
period of our observations, and in particular during the 
period of the EGRET pointing. Since electrons radiating 
optical synchrotron photons would emit y-rays by the in¬ 
verse Compton process, one inevitably predicts that the 
source had to be bright also in the EGRET band. The 
exact amount of the expected y-ray emission is of course 
model-dependent (this is the reason why simultaneous ob¬ 
servations are so important), but we can restrict the ex¬ 
isting possibilities by noting that the absence of emission 
lines in the S5 0716+714 spectrum makes inverse Compton 
off photons produced externally to the jet less likely. Eor 
S5 0716+714, a simple homogeneous and one-zone model 
is appropriate to describe the optical/y-ray behaviour, 
because the regions emitting both these frequency bands 
are cospatial. This remains true in more realistic inhomo¬ 
geneous models, which can better explain the entire elec¬ 
tromagnetic output of the source, including the radio, far 
infrared, and X-ray bands. These models can be put on 
test by simultaneous monitoring in the above mentioned 
bands. In fact, a key ingredient of the inhomogenous jet 
model we used is that smaller jet regions, located at short 


distances from the jet apex, are responsible for the high 
energy steep tail of the X-ray synchrotron spectrum, while 
larger and outer regions mainly contribute to the radio- 
to-IR band by synchrotron, and to hard X-rays by self- 
Compton emission. Perturbances travelling along the jet 
could be mapped in the frequency space, with shorter vari¬ 
ability timescales corresponding to higher synchrotron and 
inverse Compton frequencies. In addition, there should be 
time delays in the light curves at different frequencies. 

If the seed photons to be upscattered in energy 
are local synchrotron photons, then we predict that S5 
0716+714 should have reached, during the February 1995 
EGRET observation, its brightest ever detected y-ray 
state, with a flux a factor 3-10 higher than in the 1992- 
1993 EGRET observing periods. Furthermore, the simple 
SSC models presented here assume that the bulk of the 
optical and y-ray luminosities are produced in the same 
region, and therefore the flux in both bands should vary 
together, without delay. In particular, we observed a de¬ 
cline of the optical flux on a time scale of a few days, and 
we expect that a similar trend should also be found in the 
y-ray band if the statistical quality of the data is good 
enough. If the y-ray emission is not linked with the de¬ 
cline of the bulk of the optical emission, one is still left 
with the interesting alternative that the y-rays are con¬ 
nected with the optical fast flares. On the other hand, the 
absence of any optical-y-ray correlation would be a severe 
problem for all simple SSC models. 

References 

Antonucci, R.R.J., Hickson, P., Olszewski, E.W. & Miller, J.S., 
1986, AJ 92, 1. 

Barlow, 1989, Statistics, J. Wiley & Sons, Chapter 8 
Biermann, P.L., et al. 1981, ApJ 247, L53 
Biermann, P.L., Schaaf, R., Pietsch, W., Schmutzler, T., 
Witzel, A. & Kiihr, H., 1992, A&AS 96, 339. 



Ghisellini et al.: Optical-IUE observations of S5 0716+714... 


13 


Blandford, R.D., 1990, in Active Galactic Nuclei, eds. T.J.L. 

Gourvoisier, & M. Mayor (Berlin) 

Blandford, R.D. & Levinson, A., 1995, ApJ 441, 79 
Blandford, R.D., 1993, in Gompton Gamma Ray observatory, 
ed. N. Gehrels & M. Friedlander (New York: American In¬ 
stitute of Physics), 533 

Bohlin, R., Harris, A.W., Holm, A.V., & Gry, C., 1990, ApJS 
73, 413 

Gappi, M., Gomastri, A., Molendi, S., Palumbo, G.G.G., Della 
Geca, R. & Maccacaro, T., 1994, MNRAS 271, 438 
Gardelli, J.A., Glayton, G.G., Mathis, J.S. 1989, ApJ 345, 245. 
Gassatella, A., Gonzalez-Riestra, R., Imhoff, G., Oliversen, N., 
& Lloyd, G., 1992, A&A 256, 309 
Gelotti, A., Maraschi, L. & Treves, A., 1991, ApJ, 377, 403 
Giliegi, P., Bassani, L. & Garoli, E., 1993, ApJSS 85, 111 
Gomastri, A., Molendi, S. & Ghisellini, G., 1995, MNRAS 277, 
297 

Dermer, G. & Schlickeiser, R., 1993, ApJ 416, 458 
Dondi, L. & Ghisellini, G., 1995, MNRAS 273, 583 
Eckart, A., Witzel, A., Biermann, P., Johnston, K.J., Simon, 
R., Schalinski, G. & Kiihr, H., 1986, A&A 168, 17 
Edelson, R., Pike, G. F., Saken, J. M., Kinney, A., & Shull, J. 

M., 1992 ApJS, 83, 1 
Fichtel, G.E. 1994, ApJS 90, 917 

Falomo, R., Treves, A., Ghiappetti, L., Maraschi, L., Pian, E., 
& Tanzi, E. G. 1993, ApJ 402, 532 
Ghisellini, G., 1989, MNRAS 238, 449 

Ghisellini, G., Maraschi, L. & Treves, A., 1985, A&A 146, 204 
Ghisellini, G. & Maraschi, L., 1989, ApJ 340, 181 
Ghisellini, G, Padovani, P, Gelotti, A. & Maraschi, L., 1993, 
ApJ 407, 65 

Ghisellini, G. & Madau, P., 1996, MNRAS 280, 67 
Hartman, R.G., Webb, J.R., Marscher, A.P., et al., 1996, ApJ 
461, 698 

Heeschen, D.S., Krichbaum, T.P., Schalinski, C.J. & Witzel, 
A., 1987, AJ 94, 1493 

Heidt, J. & Wagner, S.J., 1996, A&A 305, 42 
Horne, K., 1986, PASP 98, 609 
Impey, G.D. & Neugebauer, G., 1988, AJ 95, 307 
Kiihr, H., Pauliny-Toth, I.I.K., Witzel, A. & Schmidt, J., 1981, 
AJ 86 , 854 

Kiihr, H., Witzel, A., Pauliny-Toth, I.I.K. & Nauber, U., 1981, 
A&AS 45, 367 

Landolt, A.U., 1992, AJ 104, 340 

Lin, Y.G., Bertsch, D.L., Dingus, B.L., et ah, 1995, ApJ 442, 
96 

Maraschi, L, Ghisellini, G. & Gelotti, A., 1992, ApJ 397, L5 
Maraschi, L. & Rovetti, F., 1995, ApJ 436, 79 
Massaro, E., Nesci, R., Maesano, M., et al., 1996, A&A 314, 
87 

Massaro, E. & Trevese, D., 1996, A&A 312, 810 
McBreen, B., 1979, A& A 71, L19 

Perley, R.A., Fomalont, E.B. & Johnston, K.J., 1980, AJ 85, 
649 

Perley, R.A., Fomalont, E.B. & Johnston, K.J., 1982, ApJ 255, 
L93 

Pian, E. & Treves, A., 1993, ApJ 416, 130 
Pian, E., Urry, M.G., Treves, A., et al., 1997, ApJ, in press. 
Quirrenbach, A., Witzel, A., Krichbaum, T., Hummel, C.A., 
Alberdi, A. & Shalinski, G., 1989, Nat 337, 442 


Quirrenbach, A., Witzel, A., Wagner, S., et ah, 1991, ApJ 372, 
L71 

Quirrenbach, A., Witzel, A., Krichbaum, T.P., et al., 1992, 
A&A 258, 279 

Reich, W., Steppe, H., Schlickeiser, R., Reich, P., Pohl, M., 
Reuter, H.P., Kanbach, G & Schonfelder, V., 1993, A&A 
273, 65 

Riecke, G.H. & Lebofski, M.J., 1985, ApJ 288, 618 
Schalinski, G.J., Witzel, A., Kircbaum, T.P., Hummel, G.A., 
Quirrenbach, A. & Johnstone, K.J., 1992, in Variability in 
Blazars, ed. E. Valtaoja & M. Valtonen, (Gambridge Univ. 
Perss), p.221 

Shull, J.M. & Van Steenberg, M.E., 1985, ApJ 294, 599 
Sikora, M., Begelman, & Rees, M.J., 1994, ApJ 421, 153 
Sillanpaa, A., Takalo, L.O., Lehto, H.J., et al., 1996, A& A, 
305, L17 

Steppe, H., Liechti, S., Mauersberg, R., Kompe, C., Brunswig, 
W. & Ruiz-Moreno, M., 1992, A&AS 96, 441 
Steppe, H., Paubert, G., Sievers, A., et ah, 1993, A&AS 102, 
611 

Stickel, M., Fried, J.W. & Kiihr, H., 1993, A&AS 98, 393 
Stickel, M., Padovani, R, Urry, G.M., Fried, J.W. & Kiihr, H., 
1991, ApJ 374, 431 
Stone, R.P.S., 1977, ApJ 218, 767 

Takalo, L.O., Sillanpaa & Nilsson, K., 1994, A&AS 107, 497 
Takalo, L.O., Sillanpaa, A. Letho, H.J., et ah, 1996, A& ASS, 
120, 313 

Thompson, D.J., Bertsch, D.L., Esposito, J.A., et ah, 1995, 
ApJS 101, 259 

Tosti, G., Pascolini, & Fiorucci, M., 1996, PASP 108, 706 
Urry, G.M. & Padovani, P., 1995, PASP 107, 803 
Urry, G.M., & Reichert, G. 1988, lUE NASA Newsletter, 34, 
96 

von Montigny, G., Bertsch, D.L., Ghiang, J., et ah, 1995, ApJ 
440, 525 

Wagner, S.J., Sanchez-Pons, F., Quirrenbach, A. & Witzel, A., 
1990, A&A 235, LI 

Wagner, S.J., 1992, in “X-ray emission form AGN and GRB”, 
MPE proceedings, p. 97 

Wagner, S.J., Witzel, A., Krichbaum, T.P. et ah, 1993, A&A 
271, 344 

Wagner, S.J. & Witzel, A., 1995, ARA&A 33, 163 
Wagner, S.J., Gamenzind, M., Dreissigacker, O., et ah, 1995, 
A&A 298, 688 

Wagner, S.J., Witzel, A., Heidt, J., et ah, 1996, AJ 111, 2187 


This article was processed by the author using Springer-Verlag 
DTgX A&A style file L-AA version 3. 



